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PSR B1259-63/LS 2883 is a binary system consisting of a 48 ms pulsar orbitting around a Be star with an orbital 
period of ~ 3.4 years. The system was detected at very high energies (VHE; E > 100 GeV) by the High Energy 
Stereoscopic System (H.E.S.S.) during its periastron passages in 2004 and 2007. Here we present new H.E.S.S. 
observations corresponding to its last periastron passage, which occurred on December 15th 2010. These new 
observations partially overlap with the beginning of a spectacular gamma-ray flare reported by the Fermi-LAT. 
The H.E.S.S. observations show both flux and spectral properties similar to those reported in previous periastron 
passages, without any signature of the emission enhancement seen at GeV energies. A careful statistical study 
based on the Fermi and H.E.S.S. lightcurves leads to the conclusion that the GeV and TeV emission during the 
flare have a different physical origin. This conlusion, in turn, allows to use Fermi-LAT measurements of the 
GeV flux as upper limits for the modeling of the VHE emission. 



1. Introduction 



The class of very high energy (VHE; E > 100 GeV) 
7-ray binaries comprises only a handful of known ob- 
jects in our Ga laxy: LS 5039 [Aharonian et al.|2005a] , 
LS I +61 303 |Albert et al.||2006| , PSR B1259-63/LS 



2883 Aharonian et al. 2005b and HESS J0632+057 



Aharonian et al.|2007 , the first binary primarily dis 



covered at TeV energies. Only for the source presented 
in this paper, PSR B1259-63/LS 2883, the compact 
companion is clearly identified as a pulsar, making it a 
unique object for the study of the interaction between 
pulsar and stellar winds and the emission mechanisms 
in such systems. 

It consists of a rapidly rotating pulsar with a spin 
period of ~ 48 ms and a spin-down luminosity of 
~ 8 x 10 35 erg/s in a highly eccentric (e = 0.87) orbi t 
around a massive Be star Johnston et al. 1992a|b 



The pulsar moves around the companion with a pe- 
riod P or b = 3.4 years (1237 days). The luminosity of 
the companion star LS 2883 is = 2.3 x 10 38 ergs -1 

Because of its fast rota- 



Negueruela et al. 



2011 



tion, the star is oblate with an equatorial radius of 
R eq = 9.7 Rq and a polar radius of i? po io = 8.1 Rq. 
This leads to a strong gradient of the surface temper- 
ature from T eq « 27, 500K at the equator to T polc ss 
34, 000A" at the poles. The optical observations also 
suggest that the system is located at the same distance 
as the star association C en OBI at d = (2.3 ±0.4) kpc 
Negueruela et al.]|2011 . The companion Be star fea- 



tures an equatorial disk that is believed to be inclined 
with respect to the pulsar's orbital plane Johnston 
et al.||1992al |Negueruela et al.||201l| |Melatos et~aL 



1995 in a way that the pulsar crosses the disk twice 
in each orbit just before (~ 20 days) and just after 
(~ 20 days) the periastron. 

Since its discovery in 1992, PSR B1259-63/LS 2883 
is constantly monitored by various instruments at all 
energy bands. The properties of the radio emission 
differ depending on the distance between the pulsar 



and the star. Radio observations Johnston et al. 1999 



Connors et al. 2002 



Johnston et al. 2005 show that 



when the pulsar is far from the periastron, the ob- 
served radio emission consists only of the pulsed com- 
ponent, whose intensity is almost independent of the 
orbital position. But closer to the periastron, starting 
at about t p — 100 d, where t p is the time of perias- 
tron, the intensity starts to decrease up to the com- 
plete disappearance approximately at t p — 20 d. This 
is followed by an eclipse of the pulsed emission for 
about 35-40 days when the pulsar is behind the disk. 
In contrast, a transient unpulsed component appears 
and sharply rises to a level more than ten times higher 
than the flux density of the pulsed emission far from 
the periastron. After the disk crossing the unpulsed 
emission shows a slight decrease with another increase 
around t p + 20 d at the second crossing of the disk. 

PSR B1259-63/LS 2883 is very well covered by X- 
ray observations carried out with various instruments 
such as ROSAT |Cominsky et al.|1994| , ASCA [Kaspi 
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et al.||1995| [Hirayama et al.|1999| , INT EGRAL |Shaw 



ct al. 2004 



etaU 2009 



Suzaku, Swift, Chandra 
and XMM-Newton 



Chernyakova 



Chernyakova et al. 



2009 2006 . X-ray observations did not show any X 
ray pulsed emission from the pulsar. Unpulsed non- 
thermal radiation from the source varies with orbital 
phase in both flux and spectral index. Similarly to 
radio measurements, the enhancement of the flux oc- 
curs shortly before and shortly after the periastron. 
Unambiguously, the enhancement of the non-thermal 
emission results from the interactions of the pulsar 
wind with the circumstellar disk close to the perias- 
tron passage. 

PSR B1259-63/LS 2883 was observed by H.E.S.S 

2004 

and 2007 



around the periastron passag es in 
et al.||2005b 



Aharonian et al. 2009 , lead 



Aharonian 



ing to a firm detection on both occasions. In 2004, 
PSR B1259-63/LS 2883 was observed mostly after the 
periastron, in 2007 mostly before it. Therefore, the 
repetitive behaviour of the source, i.e. the recurrent 
appearance of the source near periastron at the same 
orbital phase, with the same flux level and spectral 
shape of the emission, was not precisely confirmed, 
since the observations covered different orbital phases. 
However, the similar dependence of the flux on the 
separation distance between the pulsar and the star 
for both periastron passages provides a strong indica- 



tion for the repetitive behaviour Kerschhaggl 2011 



PSR B1259-63/LS 2883 was not detected in observa- 
tions performed far from periastron in 2005 and 2006, 
which comprised 8.9 h and 7.5 h of exposure, respec- 
tively. 

Observations around the 2004 and 2007 periastron 
passages showed a variable behaviour of the source 
flux with time. A combined light curve of those two 
periastron passages indicates two asymmetrical peaks 
around periastron with a significant decrease of the 
flux at the periastron itself. Peaks of the TeV emission 
roughly coincide with the flux enhancement observed 
in other wavebands as well as with the eclipse of the 
pulsed radio emission, which indicates the position of 
the circumstellar disk. This coincidence suggests that 
the TeV emission from PSR B1259-63/LS 2883 may 
be connected to the interaction of the pulsar with the 
disk. 



The pulsed radio emission was monitored with Parkes 
telescope revealing an eclipse of the pulsed signal last- 



ing from t p 



16 d to i p +15 d. 



Radio emission from 



PSR B1259-63 at frequencies between 1.1 and 10 GHz 
was observed using the ATCA array in the period 
from t p — 31 d to t p + 55 d. The detected unpulsed 
emission around the periastron passage showed a be- 
haviour similar to the one observed during previous 
observations. The X-ray energy band was covered by 
three instruments: Swift, Suzaku and XMM-Newton. 
Observations confirmed the 1-10 keV lightcurve shape 
obtained in previous periastron observations, showing 
a rapid X-ray brightening starting at about t p — 25 
d with a subsequent decrease closer to periastron and 
a second increase of the X-ray flux after periastron 
[Abdo et aLpjllj 



Observations of the binary system PSR B1259- 
63/LS 2883 at HEs were performed using the Large 
Area Telescope (LAT) on board of Fermi. The data 
taken around the periastron passage were analysed by 



two independent working groups Abdo et al. 2011 



Tarn et al. 2011 , yielding similar results for the flar- 
ing period (see below) , although there are some dis- 
crepancies related to the first detection period close to 
the periastron passage. Those differences do not af- 
fect however the conclusions drawn in this paper. The 
source was detected close to periastron with a very low 
photon flux above 100 MeV of about (1 - 2) x 10~ 7 
cni _2 s _:l . After the initial detection, the flux de- 
creased and the source was below the detection thresh- 
old until 14th of January, t p + 30 d, when a spectacu- 
lar flare was detected which lasted for about 7 weeks 
with an average flux ~ 10 times higher than the flux 
detected close to the periastron 



Abdo et al. 2011 



The highest day-averaged flux during the flare almost 
reached the estimate of the spin-down luminosity of 
the pulsar. The HE emission around the periastron as 
a function of time significantly differs from the two- 
peak lightcurves observed in other wavebands. The 
flare is not coincident with the post-periastron peak 
in radio, X-rays, and VHE 7-rays. It is also much 
brighter in comparison to the first detection of GcV 
emission close to the periastron passage Abdo et al. 



2011 Tarn et al. 2011 



2. Multi wavelength Observations around 
the 2010 Periastron Passage 

H.E.S.S. observations around the most recent peri- 
astron passage which took place on 15th of Decem- 
ber 2010 were performed as part of an extended mul- 
tiwavelength (MWL) campaign including also radio, 
optical, X-ray and high energy (HE; E > 100 MeV) 
observations. A joint MWL paper is in preparation. 

The first results of the radio and X-ray monitor 
ing of the source can be found in Abdo et al. 2011 



3. H.E.S.S. Observations around 
2010/2011 Periastron Passage 

3.1. H.E.S.S. Observations Results 

PSR B1259-63/LS 2883 observations were taken in 
five nights, namely January 9/10, 10/11 , 13/14, 14/15 
and 15/16, resulting in 6.2 h of data [H.E.S.S. Cof 



laboration et al. 2013 . The source was detected 



at an 11.5 a level. A spectral analysis of the de- 
tected excess events shows that the differential en- 
ergy spectrum of photons is consistent with a sim- 
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Figure 1: Overall differential energy spectrum of the 
VHE 7-ray emission from PSR B1259-63/LS 2883 for the 
whole observation period from 9 th to 16 th of January 
2011. The solid line denotes the spectral fit with a simple 
power law. The green band represents the 1 a confidence 
interval. Points are derived for the minimum significance 
of 1.5cr per bin. Points' error bars represent la errors. 
The fi gure is taken from |H.E.S.S. Collabor ation et ah] 
|2013l. 



pie power law dN/dE = N {E/lTcV)~ r with a flux 
normalisation at 1 TeV of N = (1.95 ± 0.32 sta t ± 
0.39 syst ) x 1CT 12 TeV _1 cm~ 2 s _1 and a spectral in- 
dex T = 2.92 ± 0.30 sta t ± 0.20 syst (see Fig. [T]) with 
a fit probability of 0.64. The integral flux above 1 
TeV averaged over the whole observation period is 
F(E > I TeV) = (1.01 ± 0.18 stat ± 0.20 svs ) x 10" 12 
|H.E.S.S. Collaboration et~aI1|2013 . 



In Fig. [2j the integrated photon flux above 1 TeV 
as a function of time with respect to periastron (in- 
dicated by the dashed vertical line) is shown. The 
light curve compiles the data from all three perias- 
tron observation campaigns spanning from 100 days 
before to 100 days after the periastron. The ob- 
served flux from the 2010/2011 observation campaign 
(in green on Fig. [2]) is compatible with the flux de- 
tected in 2004 at the similar orbital phases. Although 
the 2011 observations do not exactly overlap with the 
orbital phases of previous studies, they cover the gap 
in the 2004 data post-periastron light curve and the 
integrated flux follows the shape of the light curve, 
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Figure 2: Integrated photon flux above 1 TeV as a 
function of the time with respect to the periastron 
passage indicated with the vertical dashed line. The 
corresponding orbital phases (mean anomaly) are shown 
on the upper horizontal axis. The data from the 2004 
(blue squares) Aharonian et al. 2005b , 2007 (red 



triangles) |Aharonian et al.||2009| , and 2011 (green 
circles) observation campaigns are shown. For the 2004 
and 2011 data the flux is shown in daily bins while for 
the 2007 data the flux is shown in monthly bins for 
clarity. The figure is taken from |H.E.S.S. C ollaborat ion | 
[eFldTl f20T3l . 



yielding a stronger evidence for the repetitive be- 
haviour of the source. Moreover, the spectral shape 
of the VHE 7-ray emission from PSR B1259-63/LS 
2883 around the 2010/2011 periastron passage is sim- 
ilar to what was observed during previous periastron 
passages. The photon index of 2.92 ± 0.25 sta t ± 0.2 syst 
inferred from the 2011 data is well compatible with 
previous results H.E.S.S. Collaboration et al. 112013 . 



3.2. Search for the GeV Flare 
Counterpart at VHEs 

The absence of the flux enhancement during the 
GeV flare at radio and X-ray wavebands indicates that 
the GeV flare may be created by physical processes 
different from those responsible for the emission at 
other wavelengths. If one assumes that HE and VHE 
emission are created according to the same scenario, 
i.e. the same acceleration and radiation processes and 
sites, then a flux enhancement of the same magni- 
tude as observed at HEs should be also seen at VHEs. 
The VHE post-periastron data obtained with H.E.S.S. 
around the 2004 periastron passage do not show any 
evidence of a flux outburst at orbital phases at which 
the GeV flare is observed. However, the H.E.S.S. ob- 
servations around the 2004 periastron passage do not 
comprise the orbital phase when the GeV flare starts. 
Moreover, to compare H.E.S.S. 2004 data with the 
GeV flare observed after the 2010 periastron passage, 
one has to assume that the GeV flare is a periodic 
phenomenon, which may not be the case. 
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Table I Analysis results of the H.E.S.S. data for the pre-flare and flare periods. 



Dataset 


Livetime [h] 


Significance [a] 


1 




AT Tin — 12 rpT T- 

JMo [10 lev 


cm s ] 


t lux(£/ > 1 lev J 
















[10- 12 cm^s- 1 ] 


Pre-flare 


2.65 


7.4 


2.94 ± 0.52 stat 


± 0.20 syst 


2.15±0.56 st at 


± 0.43 syst 


1.11 ± 0.29 sta t ±0.22 syst 


Flare 


3.59 


8.5 


3.26 ± 0.49 stat 


± 0.20 syst 


1.81 ± 0.39 stat 


± 0.36 S yst 


0.80 ± 0.22 stat ±0.16 syst 



_, 3r 



and a x 2 -to-NDF ratio o f 0.73/1, which corresponds to 



,i. H.E.S.S. 




31 32 
Days from periastron 



Figure 3: Integrated photon fluxes above 1 TeV for the 
pre-flare and flare periods (see text) are shown as black 
filled boxes. The dashed horizontal line shows the best fit 
with a constant. The HE data points above 0.1 GeV as 

are shown as red filled 



2011 



reported by Abdo et al. 
circles. The flare start date is indicated by the dashed 
vertical line. The left axis indicates the units for the 
VHE flux and the right (red) axis denotes the units for 
the HE flux. The figure is taken from 



H.E.S.S. Collaboration et al. 2013 



The H.E.S.S. data taken between 9th and 16th 
of January in 2011 provide a three-day overlap in 
time with the GeV flare. Therefore, it is possible 
to directly study any flux enhancement in the VHE 
band on the time scale of the HE flare. To im- 
prove the sensitivity of the variability search the whole 
period of the H.E.S.S observations was divided into 
two time periods of almost equal length: "pre-flare" 
(i p + 26 d to t p + 29 d; before the HE flare) and "flare" 
(i p + 30 d to t p + 32 d; during the HE flare) . These two 
datasets were analysed independently and revealed 
similar fluxes and significance levels (see Table [I]) . A 
spectral analysis of the two datasets shows that both 
spectra are consistent with a simple power law, yield- 
ing similar values of the spectral index (see Table [I]) . 
The two spectral indices are consistent with the one 
obtained for the total dataset 
et^[2013 



H.E.S.S. Collaboration! 



A fit of the flux as a function of time with a con- 
stant showed no indication of the flux variability. The 
fit resulted in a mean flux of (0.91 ± 0.18) x 10~ 12 
enr 2 s _1 (black horizontal dashed line in Fig. [3] top) 



2013 



a x probability of 0.39 H.E.S.S. Collaboration et al. 



To firmly reject the hypothesis of the GeV flare 
counterpart at VHEs, i.e. the flux enhancement of 
the same magnitude as GeV flare, the flare coefficient 
k was introduced as the ratio of the fluxes during the 
flare period and the pre-flare period H.E.S.S. Collab 



oration et al]|2013| . The ratio of the HE (E > 0.1 



GeV) flux averaged over the three-day interval be- 
tween (t p + 30 d) and (t p + 32 d) to the upper limit 
on the HE pre-flare emission (see Fig. [3] (Top)) 
yields a lower limit on the HE emission flare coeffi- 
cient khe > 9.2. The 99.7 % confidence level (equiv- 
alent of 3cr) upper limit on the VHE flare coefficient 
was estimated using the profile likelihood method (see 
H.E.S.S. Collaboration et al.| [2013] for details). The 



obtained upper limit of Kgg.7% < 3.5 is lower than the 
observed lower limit on «he- This result suggests that 
the HE flare emission has a different nature than the 
VHE emission. 



4. Modeling of the VHE Emission 

Several models have been proposed to explain the 
VHE emission from the source. In a hadronic sce- 
nario, the VHE 7-ray emission could be produced by 
the interaction of the ultrarelativistic pulsar wind par- 
ticles with the dense equatorial disk outflow with sub- 
sequ ent production of 7r° pions and hence VHE 7- 
rays Kawachi et al. 112004 Neronov and Chernyakova 



2007 



However, the detection of the source before 
the expected disk passage in 2007 casts doubts on 
the hadronic scenario, suggesting that the VHE emis- 
sion sh ould be created at least partly by leptonic pro- 
cesses lAharonian et al. 2009]. Within the leptonic 



scenario, VHE emission from PSR B1259-63/LS 2883 
is explained by the inverse Compton (IC) scatte ring 
of shock-accelerated electrons on stellar photo ns |Ta-| 
vani and Arons||1997[ |Kirk et al.||1999[ |Dubus||2PjOl3aT 



Khangulyan et al.||2007 



In this paper for the modeling of the VHE emission 
the leptonic scenario is considered. Since the flaring 
HE emission was produced in a different scenario than 
the VHE emission, the HE flux measured by Fermi- 
LAT can be used to constrain parameters which de- 
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scribe the electron population. Few possible expla- 
nations for the nature of the HE flare are discussed 



in the literature Khangulyan et al. 2012 


Bogovalov 


et al.||2008 


Dubus et al.||2010 Kong et al. 


2012 , but 



their discussion is beyond the scope of this paper. 

Since at VHE the system was observed for a short 
period and the lightcurve does not show any variabil- 
ity on this time scale, the time dependence can be 
neglected, and thus, one population of electrons can 
be assumed. The energy distribution of the electron 
density is assumed to follow a power-law with an ex- 
ponential cut-off 



dN c 
e?7 



KfSy p e th 



(1) 



the 



where p is the electron spectral index, K c is 
normalisation parameter with 7 = 1, and -E max = 
wic 2 7 max is the electron cut-off energy. It is further 
assumed that this distribution of electrons already ac- 
counts for all kinds of energy losses and that these 
losses are not important on the considered time scale. 
The absorption due to the pair production of the in- 
dividual 7-rays is not taken into account in the calcu- 
lation of the resulting 7-ray flux, as it does n ot play 
an important role in the post-periastron phase Dubus 



2006b . The electrons are assumed to be distributed 
isotropically at the spherical termination shock. The 
stellar radiation at the location of the pulsar wind ter- 
mination shock plays the role of target photons. The 
target photon field has then a Planckian distribution 
with a temperature = 30000-ftT and an energy den- 
sity of 1 erg cm -3 , which corresponds to the distance 
between the star and the shock at the observed orbital 
phase, assuming that the distance from the pulsar to 
the shock is much smaller than the distance from the 
star to the shock. Although the model is rather simple 
and does not account for a number of different pro- 
cesses which take place in the system in the vicinity of 
the emission region, it can still give an understanding 
of the population of electrons required to produce the 
observed VHE 7-ray flux. 

To constrain the IC flux the Fermi upper limit 
value at 100 GeV is used. If one assumes that -Eul = 
100 GeV is far from the cut-off energy, the IC emission 
at this energy would depend only on two parameters, 
K e and p. Therefore, for each assumed p, one can 
calculate the upper limit on K e . Then, for each pair 
of p and K e one can estimate -E max by fitting the IC 
emission to the H.E.S.S. data (Fig. [4}. In Table |l!) 
upper limits on K c and corresponding £? max obtained 
from the fit of the H.E.S.S. data are shown for differ- 
ent values of the spectral index. The fit probability is 
also given in the table and is > 0.5 for all studied spec- 
tral indices. Based on the obtained parameters of the 
electron distribution an upper limit on the total en- 
ergy in electrons Wtot which is needed to produce the 
VHE emission not exceeding the HE emission upper 
limit can be calculated. It is calculated assuming the 




Figure 4: Inverse Compton emission models for different 
spectral indices of the electron energy distribution. 

H. E.S.S. data are shown with black squares and 
Fermi-LAT data are shown with red circles. 

minimum energy of the electron -Ec.min = 1 GeV. De- 
pending on the assumed spectral index it varies from 
4.0 x 10 43 erg for p = 1.7 to 1.7 x 10 46 erg for p = 2.9. 
Assuming that the spin-down luminosity of the pul- 
sar, which is about 8 x 10 35 erg/s, is transferred to 
the electron acceleration with the 100% efficiency a 
lower limit on time required to accumulate Wtot can 
be derived. The lower limit on time ranges between 

I. 5 y for p — 1.7 to 674 y for p = 2.9. For the range 
of spectral indices from 1.7 to 1.9 the lower limit on 
time is less than 3.4 y, i.e. the orbital period of the 
pulsar. This would mean that electrons accelerated at 
the shock during one orbital period may lose all their 
energy producing IC emission close to periastron and 
this process may repeat every 3.4 y. A possible exten- 
sion of the electron spectrum down to energies lower 
than 1 GeV would not change results significantly for 
p < 2.0. 

Table II Fit parameters for fixed values of the spectral 
index p 



p 


(4ttD2) ' Cm 


£max, TcV Fit prob. 


Wtot, erg 


1.7 


2.2 x 10 2 


6.2 


0.72 


4.0 x 10 43 


1.9 


5.0 x 10 3 


7.1 


0.70 


6.9 x 10 43 


2.1 


1.1 x 10 5 


8.4 


0.68 


1.5 x 10 44 


2.3 


2.6 x 10 6 


10.1 


0.65 


4.0 x 10 44 


2.5 


5.8 x 10 7 


12.8 


0.61 


1.3 x 10 45 


2.7 


1.3 x 10 9 


17.3 


0.56 


4.5 x 10 45 


2.9 


2.9 x 10 10 


26.1 


0.50 


1.7 x 10 46 



5. Summary 

The binary system PSR B1259-63/LS 2883 was 
monitored by H.E.S.S. around the periastron passage 
on 15th of December 2010. The observed flux and 
spectral shape agree well with what was measured 
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during previous periastron passages. The observations 
were performed at similar orbital phases as around the 
2004 periastron passage, strengthening the evidence 
for the repetitive behaviour of the source at VHEs. 

H.E.S.S. observations were part of a joint MWL 
campaign that also included radio, optical, X-ray, and 
HE observations. A spectacular flare observed at HEs 
with Fermi LAT overlapped in time with the H.E.S.S. 
observations. A careful statistical study showed that 
the HE flare does not have a counterpart at VHEs, in- 
dicating that the HE and VHE emissions are produced 
in different physical scenarios. The measured HE 
flux is then used to constrain the emission produced 
through the inverse Compton scattering of shocked 
relativistic electrons of stellar photons which is re- 
sponsible for the VHE flux. It is shown that for the 
range of assumed electron spectral indices from 1.7 
to 1.9 the lower limit on time is less than the orbital 
period of the pulsar. 
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